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1. Introduction 

The flavour physics program plays a dominant role in testing the Standard Model (SM) and 
searching for New Physics (NP), providing information which is complementary to direct searches 
in colliders. NP effects could be unveiled through the observation of deviations from the SM via 
high-precision measurements of low-energy observables in high-luminosity experiments. Indeed, 
there are several measurements of flavour observables for which there is a 2 — 3a difference from 
SM predictions. These include sin(2j3), the like-sign dimuon charge asymmetry, B leptonic decays, 
unitarity triangle (UT) fits, and, more recently, the ratios of branching fractions of B semileptonic 
decays to a D and a z over the corresponding decays to a / = e,pL. 

In order to test the Cabibbo-Kobayashi-Maskawa (CKM) description of the experimentally 
measured CP-violating and flavour changing processes, and investigate the origin of the discrep- 
ancies mentioned above, we need a determination of the weak matrix elements involved in those 
processes with matching precision. In many cases, lattice QCD can provide those non-perturbative 
theoretical inputs from first principles and with errors at a few per-cent level. Accuracy in lattice 
calculations requires control over all the sources of systematic error. In particular, it is essential 
to take into account vacuum polarization effects in a realistic way, i.e., including up, down and 
strange sea quarks on the gauge configurations' generation. The up and down quarks are usually 
taken to be degenerate, so those simulations are referred to as Nf = 2 + 1. Two lattice collabora- 
tions (FNAL/MILC and ETMC) are now generating configurations which also include the effects 
of charm quarks on the sea, Nf = 2 + 1 + 1. The first preliminary results for flavour quantities on 
those configurations are starting to appear. 

In the next Sections I will discuss the latest results for non-perturbative quantities relevant 
for flavour studies from lattice QCD calculations with all sources of systematic error addressed. 
Among other things, that means that I will focus on simulations with Nf = 2 + 1 sea quarks. Due to 
the increasing number of results available for some quantities from different collaborations working 
with Nf = 2 + 1 configurations, it is important to perform averages of those results that can be used 
by non-experts. There are already two groups working on initiatives to compile results and provide 
averages: FLAG [0] and Laiho-Lunghi-Van de Water [§]. A more ambitious project including 
members of those two averaging collaborations and new collaborators from Europe, Japan, and 
the US, the Flavour Lattice Averaging Group, FLAG-2, expects to provide its first set of averages 
soon [J3|] . 

2. Light quarks matrix elements 

2.1 \V US \ from K leptonic and semileptonic decays 

The lattice calculation of pseudoscalar decay constants, together with experimental measure- 
ments of pseudoscalar leptonic decay widths, can be used to extract the value of the CKM matrix 
elements involved in those processes. The decay constants are easy to calculate on the lattice with 
high precision. The accuracy achieved is even higher for ratios of decay constants, for example 
fulfil, since many systematic uncertainties and statistical fluctuations cancel partially or com- 
pletely between numerator and denominator, fx/ fn has been widely studied on the lattice and the 
precision has reached the subpercent level [ffl]. The average of the published results with Nf = 2 + 1 
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is [fj] fa/fx = 1-1936 ±0.0053. Agreement between different collaborations also provides a good 
check of lattice methodologies. The phenomenological interest of this quantity stems from the fact 
that it is related to the ratio |V us |/|Kd| and experimental measurements of K and 71 leptonic decay 
widths, so it can be used to extract a value for \V US \. Using the experimental average in and the 
lattice average for fx/f% above, we get \V US \ = 0.2252(11). 

Precise determinations of \V US \ provide stringent tests of first-row unitarity and give informa- 
tion about the scale of NP This CKM parameter can also be extracted from experimental 
data on K semileptonic decay rates, given the value of the vector form factor at zero momentum 
transfer, f+ n (§). The parameter ff n {0) is also one of the inputs used in UT analyses. The only 
unquenched lattice calculations available for this form factor 1 are the ones by the RBC/UKQCD 
and ETMC collaborations [§]. The average of the two calculations is [g] ff n (0) = 0.9584(44), 
which gives \V US \ = 0.2257(11), in perfect agreement with the value extracted from leptonic decays 
and also in agreement with unitarity. Although these errors are already very small, there is a lot of 
room for improvement: more lattice spacings, more sophisticated extrapolation methods, physical 
light quark masses, ... The RBC/UKQCD, FNAL/MILC, and JLQCD collaborations have already 
presented preliminary results incorporating some of those improvements in [^J. 

2.2 K Q -K° mixing 

The information coming from neutral Kaon mixing provides one of the most stringent con- 
straints in UT analyses. Until a few years ago, the limiting factor to exploit that consuaint was the 
uncertainty in the calculation of the hadronic matrix element which encodes the non-perturbative 
physics of the process, parametrized by the bag parameter Bk- This is no longer true thanks to 
several Nf = 2 + 1 lattice QCD calculations, four of them in the last year JT0|], which have reduced 
the error of Bk under 2%. The average of all Nf ■ = 2 + 1 calculations is [Q] Bk = 0.7643(97). The 
dominant errors in the \ek \ constraint in UT analyses are now the uncertainty in |V C £,|, which enters 
as the fourth power in the theoretical prediction for |%|, and the error in the NNLO perturbative 



coefficient r\ cc Jl This year we have also had the first unquenched calculations of the matrix 
elements of the four-fermion operators contributing to K° — K° in extensions of the SM by the 
RBC/UKQCD Collaboration @ (Nf = 2 + 1 calculation at a single lattice spacing) and by the 
ETMC Collaborarion JT3| ] (Nf = 2 calculation in the continuum limit). 



3. Heavy quark phenomenology 

3.1 Heavy-light decay constants 

The decay constant of the D s meson, fu s , has generated a vivid interest in both experimentalists 
and lattice QCD phenomenologists in the last years, especially since 2007, when an important 
reduction in the errors of the lattice calculation of fo s made possible to observe a 3.8 a disagreement 
between the lattice and experimental averages. Later improvements in both the theoretical and the 
experimental sides, however, shift both numbers and reduced the difference. Nevertheless, the most 

'After this conference, the FNAL/MILC collaboration has provide another result [Q], compatible with those by the 
RBC/UKQCD and ETMC collaborations, but with a larger central value and smaller errors. The FNAL/MILC result 
implies |V,„| = 0.2238 ±0.0009 ±0.0005. 



3 



Flavour physics from lattice QCD 



Elvira Gamiz 



precise lattice calculations, including those in a preliminary stage, tend to give smaller values of 
fo s than experiment. In particular, two Nf = 2 + 1 + 1 preliminary results from the ETMC and the 
FNAL/MILC collaborations, are around (245-250) MeV [0]. The lattice averages, which only 
include the Nf = 2 + 1 calculations already complete, are [|[ [TJ|, |l6|] 

f D = 209.2(3.9)MeV , f Ds = 248.6(3.0)MeV . (3.1) 

The value of fo agrees with experiment [|7]|, fo = 206.7 ±8.9, while fo s is around 2a lower than 
the experimental average @], f Ds = (260.0 ± 5.4)MeV. The preliminary BESIII result in [||], 
fo = (203.9 ± 6.0)MeV, agrees well with the CLEO result and the lattice average. 

UT fits are very sensitive to fg and different processes with potential to show up NP effects 
depend on or fg s , so any improvement in the decay constants calculations as well as in the 
understanding of the \V u b\ exc ' /\V u b\ mc ' disagreement is very important. There have been three lat- 



tice Nf = 2 + 1 calculations of this parameter in the last two years by the HPQCD Q20i |21| ] and 



the FNAL/MILC [22] collaborations, which have reduced the error at the 2.5% level. The small- 



est error in plj ] is achieved by determining the ratio /bJ/b using a non-relativistic description 



(NRQCD) of the b quarks together with the fg s determination in [|20|], which employs relativistic 
b quarks. By using the ratio /bJ/b the dominant systematic errors associated with the effective 
NRQCD description are partially cancelled, so this determination is nearly free of uncertainties 
due to an effective description of the b— quark. The average values of fs and fg s from these three 
calculations are [Q] 

f B = 190.6(4.7)MeV , f Bs = 227.6(5.0)MeV . (3.2) 



The direct comparison of the results in Eq. (3.2) with experimental measurements of the B leptonic 
decay width is problematic due to the need of the value of the CKM matrix element \V u b\ (whose 
inclusive and exclusive determinations disagree at the 3a level) and the ~ 2a disagreement of 



BaBar [23] and Belle's [24] measurements. Nevertheless, Belle new result seems to alleviate the 



tension between theory and experiment previously observed. 

3.2 B semileptonic decays 

There exist ~ 2 — 3a tensions between the inclusive and the exclusive determinations of 
both \V u b\ and \V c b\- In addition to experimental measurements on, for example, B — >• nlv and 
B — > D{D*)lv, respectively, the exclusive determination of those CKM elements need as input 
form factors that can be calculated with high precision using lattice QCD techniques. There have 
not been new lattice QCD calculations of the form factors describing the B — > nlv decay since 



2008 [|25j], although several analyses are in progress [|26J|. Using the lattice results in Ref. [ |25| ] and 
the latest experimental data [27[l, |V H fc I exclusive = (3.23 ±0.30) • 10~ 3 . 



For the exclusive determination of \V c b\, the state-of-the art calculation of the relevant form 
factors is the FNAL/MILC analysis in p8|], which studies the decay B — > D*lv at zero recoil. 
The updated result is \V c b\ = (39.7 ±0.7 ±0.7) • 10~ 3 , where the first error is experimental and 
the second one the uncertainty in the lattice calculation of the form factors. The FNAL/MILC 
collaboration is doing an extensive study of both B — > D*lv and B — > Dlv decays at zero and non- 
zero recoil p9[|, providing two independent modes for the extraction of \V c b\. This study will also 
provide checks of the shape of the form factors, in addition to \V c b\- 
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The FNAL/MILC collaboration recently analyzed a subset of their B — > Dlv data to calculate 
the ratio of branching fractions R(D) = 38r(B ->■ Drv)/&r(B ->■ Dlv) = 0.316(14) Q. Their 
value is ~ 1 .7a smaller than the recent experimental measurement by the BaB ar collaboration [[H]] . 
They found that the value of the ratio is very sensitive to differences in the scalar form factor, so 
one should be cautious in using indirect estimates of the form factors to constrain NP models in 
other decay channels such as B — > D*tv [30]. Given the present tensions, not only for R(D) but 
for R(D*), and the possible indications of NP that could be extracted from a combined analysis of 
both set of decays, unquenched lattice QCD calculations of those two ratios should be a priority. 
Together with a determination of \V c b\ and the shape of the form factors describing B — > D(D*)lv 
decays, the final FNAL/MILC analysis including the complete set of data will also provide an 
improved determination of R{D) as well as a calculation of R(D*). Current experimental measure- 
ments of these ratios are statistics-limited, so Belle II should significantly reduce the errors of those 
measurements. The Cambridge group and the FNAL/MILC collaboration have also presented pre- 
liminary results for the different form factors describing the rare decays B — > K(K*)l + l~ [32, 33]. 



3.3 D semileptonic decays 

The semileptonic modes D — > Klv and D — > nlv allow us not only to extract the values of the 
CKM matrix elements |V cs ( d )|, but to test lattice QCD techniques and methodology by comparing 
the shape of the corresponding form factors, f+ K<yK \q 2 ), with experimental data. The method 
developed by the HPQCD collaboration that employs a Ward identity to relate the matrix element 
of a vector current to that of the corresponding scalar current, and the use of highly improved lattice 
discretizations to treat the charm quarks relativistically, have allowed a reduction of the errors in 
the lattice calculation of the form factors f?~* n (0) and f®^ K (0) from around a 10% to a 5% 



and a 2.5% respectively. The average of the results in [|34j g5J by the HPQCD and FNAL/MILC 
collaborations are 



|y c ,| = 0.961(ll)(24), |V cd | = 0.225(6) (10), (3.3) 

where the first error is from experiment and the second one is the lattice error in 

is compatible with unitarity (\V CS \ = 0.97345(16) and \V c d\ = 0.2252(7)), and it is competitive with 

the determination of \V c d\ from neutrino scattering, \V c d\ = 0.230(11) p6|]. 

There are several on-going projects that study the dependence of the form factors on the mo- 
mentum transfer, showing a good description of the experimentally measured shape. The final goal 
is performing a global fit, including all available experimental data on the decay rates and the lattice 
determination of the form factors at several values of the momentum transfer, to extract the value 
of the corresponding CKM matrix element in an optimal way as well as testing the consistency 
of the lattice description throughout the energy range. HPQCD has quoted a preliminary number 
extracted from a global fit to Belle, BaBar, and CLEO data, |V„| = 0.965(14) [0. 



3.4 Neutral B— meson mixing 

It has been argued that differences observed between measurements of some flavor observables 
and the corresponding Standard Model (SM) predictions may be due to beyond the SM (BSM) 
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HPQCD 


FNAL/MILC 


RBC/UKQCD 


§ 1.258(33) 


1.27(6) 


1.13(12) 


B Bs /B Bd 1.05(7) 


1.06(11) 




HPQCD: f B .ffi,= 


= 266(6)(17) MeV, 


B Bs = 1.33(6) 


HPQCD: f Bd yfi^ = 


216(9)(13) MeV, 


B Bi = 1.26(11) 



Table 1: B— meson mixing parameters. S, is defined as the ratio of the parameters in the second and third 
rows. In the case where there are two errors, the first one is statistical and the second one systematic. 



physics affecting the neutral B-meson mixing processes [|38| [39j]. Although the most recent anal- 
ysis seem to indicate that there are not large BSM contribution to neutral B-meson mixing [|40|], 
the future will bring new twists, and precise calculations of the non-perturbative inputs parametriz- 
ing the mixing in the SM and beyond are necessary for a thorough understanding of quark flavor 
physics. The current status of Nf = 2 + 1 lattice calculations of the non-perturbative quantities 
parametrizing the mass differences between the heavy and the light mass eigenstates in both the B° 
and systems, as well as the SU(3) breaking ratio is summarized in Tab. [jj The HPQCD [41] 



and FNAL/MILC [ |42| ] collaborations use the same light quark formulation, but a different descrip- 
tion for the b quarks. In the exploratory study by the RBC/UKQCD collaboration in the table heavy 
quarks are static p3|]. The average of the results in Tab. |l]for E, gives the value E, = 1.251 ±0.032. 

Beyond the SM, the mixing parameters can have contributions from AB = 2 four-fermion op- 
erators which do not contribute in the SM. The matrix elements of the five operators in the complete 
basis describing AB = 2 processes, together with the Wilson coefficients for those operators cal- 
culated in a particular BSM theory and the experimental measurements of the mixing parameters, 
can provide very useful constraints on that BSM theory. There is not a final unquenched lattice 
calculation of the matrix elements of all the operators in the AB = 2 effective hamiltonian, but 
FNAL/MILC and ETMC presented preliminary results for the complete basis in p4[ 45]. Those 



preliminary results include the matrix elements needed for the determination of the decay width 
differences, AF^^, for which there is not currently an unquenched calculation in the continuum 
limit. 



The authors of Ref. [46] suggested that the branching fractions of the rare decays B q — > }X + jX' 



(for q = s,d) could be determined from the experimental measurement of the mass difference in 
the neutral B^-meson system, AM q , and the lattice calculation of the bag parameter B Bq since 

Mr{B q ^ n + n~)/AMq = [known factors] /B q (3.4) 

In order to compare experimental measurements and the theory predictions for the decay rate 
of B° s , one should include the effects of a non- vanishing \Y S p7[|. This can be done in the SM by 



rescaling the theory prediction by 1/(1 — y s ), where y s = x Bs AT s /2 Q471]. Multiplying Eq. ( |3.4| ) by 
that factor for the B° — > jU + /i~ decay and using the HPQCD determination of the bag parameters 
B Bs = 1.33(6) and B Bd = 1.26(11) together with x Bs = 1.497(15)^, r Bd = 1.519(7)/w [||], 
and ATs = 0.1 16 ± 0.019/w- 1 [||], one gets 

mr(B s -> Ai + AOk = (3.65 ±0.20) x 10~ 9 , ,%r{B d -> M+aO = (1.04 ±0.09) x 10" 10 (3.5) 
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The direct calculation of these branching fractions have become competitive with the one in 
(3.5) J4^ ] thanks to the recent improvements in the calculation of the B-meson decay constants on 
the lattice described in Sec. IA. Including the correction factor 1/(1 —y s ) for the B s — > pL + pL~ decay 
rate, and using the same inputs as in Ref. [|^] except for fg and fg s , for which I use the averages 
described in Sec. 3.1, and Tg , which I take equal to its PDG 2012 value, Tg = \ A91{\5)ps, I get 



h(B s 



(3.64±0.23) • 10~ y , mr(B d -> = (1.07 ±0.10) • 10~ 10 . (3.6) 



The agreement between the two set of numbers in ( |3.5| ) and ( |3.6| ) is excellent, giving us confidence 
on the SM prediction for these branching fractions. This is very important since LHC is approach- 
ing the SM predictions, with the first evidence for one of these two processes recently reported by 
LHCb @. The LHCb measurement [ftj], 0§r{B s -4 jU+^i ) = (3.2+} f) ■ 10~ 9 , is consistent with 



the SM prediction in Eqs. ( |3.5[ ) and (p.6|). 

Another recent contribution of lattice QCD to the study of B s — > pl + pL~ rare decays has been 
the calculation by the FNAL/MILC collaboration of form-factors ratios between the semilep- 
tonic decays B° — >■ D + l~v and B° s — > D+/~v [51]. These ratios are a key theoretical input in 
a new strategy to determine the fragmentation fractions of neutral B decays, which are needed 
for the experimental measurement of 3Sr(B s — >• jx + ll~). The result for the ratio of form factors, 
/o(w^)//o ( m lt) = 1-046(44)(15), gives a value of the fragmentation fraction which agrees well 
with that of the Dfn /D + K~ hadronic method and with LHCb's determination via a method em- 
ploying semileptonic decays. This calculation included only a subset of the B — > Dlv data which is 
being analyzed by the FNAL/MILC collaboration for the determination of \V c b\, so we can expect 
a considerable improvement when they finish analyzing their full data set. 



4. Outlook 

Lattice QCD calculations of non-perturbative parameters relevant for flavour phenomenology 
have achieved accuracies at the per-cent level for many key quantities. The agreement between 
results from different collaborations for the same quantities also allows an important check of lattice 
methods, especially in the light and charm sectors. In the bottom sector, the first Nf = 2+1 results 
using a relativistic approach are very promising. For the next two years we expect new results from 
the FNAL/MILC, ETM, and RBC/UKQCD collaborations for decay constant, B° - B° mixing, 
and B — > nlv. One of the important steps in the improvement of lattice calculations expected 
in the next year is using simulations at the physical light quark masses. This would drastically 
reduce one of the main systematic, the one associated with the chiral extrapolation. There will 
be results with physical quark masses within a year from different collaborations: FNAL/MILC, 
BMW, RBC/UKQCD, and PAC-CS. The reduction of the dominant sources of errors will make 
necessary to include some uncertainties which were subdominant until now like isospin breaking, 
electromagnetic effects, or dynamical charm quarks. Another target for lattice calculations is going 
beyond gold-platted quantities and develop methods to calculate more demanding quantities (K — >■ 
%%, weak decay to resonances, non-local operators, ...). 
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